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Abstract 

Chemical  kinetic  modelers  make  extensive  use  of  shock  tube  ignition  data  in  the 
development  and  validation  of  combustion  reaction  mechanisms.  These  data  come  from 
measurements  using  a  range  of  diagnostics  and  a  variety  of  shock  tubes,  fuels,  and  initial 
conditions.  With  the  wide  selection  of  data  available,  it  is  useful  to  realize  that  not  all  of 
the  data  are  of  all  the  same  type  or  quality,  nor  are  all  the  data  suitable  for  simple,  direct 
comparison  with  the  predictions  of  reaction  mechanisms.  We  present  here  a  discussion 
of  some  guidelines  for  the  comparison  of  shock  tube  ignition  time  data  with  reaction 
mechanism  modeling.  Areas  discussed  include;  definitions  of  ignition  time;  ignition  time 
correlations  (with  examples  taken  from  recent  n-heptane  and  iso-octane  measurements); 
shock  tube  constant-volume  behavior;  shock  tube  diameter  and  boundary  layer  effects; 
carrier  gas  and  impurity  effects;  and  future  needs  and  challenges  in  shock  tube  research. 

Introduction 

Chemical  kinetic  modelers  have  made  extensive  use  of  shock  tube  ignition  data  in 
the  development  and  validation  of  combustion  reaction  mechanisms.  These  data  are 
available  for  many  different  fuels  and  a  range  of  reaction  conditions.  With  the  wide 
selection  of  data  that  are  available,  it  is  useful  to  realize  that  not  all  the  data  are  of  all  the 
same  type  or  quality,  nor  are  all  the  data  suitable  for  simple,  direct  comparison  with  the 
predictions  of  reaction  mechanisms.  What  is  needed  is  to  match  accurate  and  well- 
defined  experimental  shock  tube  data  with  the  predictions  of  an  appropriate 
computational  model. 

One  zeroth-order  approximation  model  (ZOAM)  regularly  applied  to  reflected 
shock  wave  ignition  time  data  is  the  following:  use  published  T5,  P5  (reflected  shock 
wave  initial  temperature  and  pressure),  and  Xi  (initial  mole  fractions  of  mixture 
components);  use  a  constant  volume  or  density  constraint  (zero  velocity,  U5=0)  combined 
with  the  assumption  of  no  heat  transfer  or  mass  diffusion  in  the  computation;  and 
compare  published  ignition  times  or  species  concentration  time  histories  to  computed 
values. 


A  higher-order,  fuller  description  of  the  shock  wave  experiment  would  include  the 
temporal  and  spatial  variations  in  the  test  gas  control  volume  caused  by  boundary  layer 
growth,  attenuation  of  the  incident  shock  wave,  and  energy  release  of  the  combustion  test 
gas  mixture.  At  some  point  in  the  progress  of  a  reflected  shock  experiment  in  a  high 
concentration  mixture,  a  transition  to  detonation  can  occur  that  would  make  the 
interpretation  of  shock  tube  data  even  more  complicated. 
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Thus  every  aspeet  of  the  ZOAM  is  open  to  interpretation  and  needs  to  be  well 
defined  to  be  useful.  Initial  state  conditions  in  the  ZOAM  need  to  be  chosen  from  the 
spatially  and  temporally  varying  values:  T5(t,x),  P5(t,x),  ¥5(1, x),  and  published  T5  and  P5 
need  to  be  defined  in  this  context,  particularly  if  averaging  or  extrapolation  is  employed, 
as  is  common.  Any  vibrational  non-equilibrium  and  relaxation  processes  that  occur, 
primarily  in  the  carrier  gases  also  affect  initial  conditions.  And  while  it  is  expected  that 
the  initial  mole  fractions  of  the  reactants  are  uniform  over  the  reflected  shock  regime 
(i.e.,  Xi(t=0,x)  =  Xi(t=0)),  reactants  that  can  rapidly  decompose  during  the  passage  of  the 
lower  temperature  and  pressure  incident  shock  may  have  some  spatial  variation.  Finally 
the  ignition  time  Xign  needs  to  be  clearly  and  properly  defined.  In  many  ignition 
experiments  the  exponential  branching  process  is  sufficiently  slow  that  the  “rapid  ignition 
process”  may  take  substantially  more  time  than  ~1%  of  the  induction  time  and  a  more 
unequivocal  ignition  time  definition  is  needed.  At  the  other  extreme,  variations  in 
ignition  time  can  occur  along  the  length  of  the  shock  tube  control  volume,  especially 
during  very  energetic  ignition  events,  and  an  ignition  time  measurement  location  needs  to 
be  defined. 

The  information  that  could  be  used  to  improve  the  ZOAM  approximation  is 
generally  not  presented  in  the  literature.  This  information  could  take  the  form  of  T5(t,x), 
P5(t,x),  V5(t,x)  (velocity)  and  Xi(t,x)  time  histories  during  the  entire  ignition  event,  or 
include  detailed  fluid  mechanical/chemical  models  that  more  accurately  predict  and 
describe  the  actual  flow  state  and  condition  throughout  the  test  gas  volume.  All  this 
information  points  to  a  movement  away  from  the  ZOAM  to  a  fuller  description  of  shock 
wave  ignition  experiments.  However,  validation  tests  of  these  fuller  descriptions  have 
rarely  if  ever  been  done  and  are  likely  to  be  facility  specific;  hence  the  ZOAM  is  likely  to 
remain  important. 

Uncertainties  in  the  predictive  ability  of  the  zeroth-order  model  can  vary  widely. 
The  zeroth-order  approximation  model  is  expected  to  be  most  accurate  under  the 
following  conditions:  short  test  times  (less  than  500  |is),  large  diameter  shock  tubes 
(greater  than  10  cm),  low  initial  fuel  and  oxidizer  concentrations  (such  that  the  pressure 
change  and  energy  release  during  ignition  are  small),  an  inert  diluent  gas  such  as  argon 
(more  uniform  shock  tube  performance  than  with  diatomic  gases  such  as  nitrogen),  and 
pre-shock  pressures.  Pi,  greater  than  some  modest  minimum  pressure  (of  order  10  torr.) 

However,  experimentally,  low  initial  fuel  concentrations  can  be  difficult  to 
produce  and  measure  accurately.  As  well,  much  of  the  previous  ignition-time  work  was 
done  on  relatively  small  diameter  shock  tubes  with  pressure  transducers  as  the  indicating 
diagnostic.  This  diagnostic  gives  useful  ignition  times  only  for  high  concentrations  (of 
order  1%  fuel  or  greater).  The  requirements  of  engine  (IC,  HCCI,  Cl  and  SI)  modelers 
also  present  a  different  set  of  problems.  Engine  modelers  would  like  ignition  time  data 
for  conditions  near  to  those  of  practical  devices.  These  conditions,  with  temperatures  as 
low  as  700  K,  pressures  as  high  as  100  atm  and  stoichiometric  mixtures  of  fuel  and  air  (1- 
2%  fuel,  20%  O2,  79%  N2)  can  result  in  ignition  times  of  several  to  tens  of  milliseconds, 
and  represent  an  extension  from  the  conditions  where  shock  tube  operation  is  best 
described  with  the  ZOAM. 
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This  is  not  to  say  that  shock  tubes  are  unsuited  for  measurements  in  these 
regimes.  Shoek  tubes  are  exeellent  devices  to  perform  these  studies  and  measure  ignition 
times.  What  is  neeessary  is  an  improved  interpretation  of  the  results  of  these 
measurements  as  well  as  improved  understanding  of  the  capabilities  of  shock  tubes  and 
more  eomplete  eharaeterization  of  (faeility-speeific)  shock  tube  behavior  and  flow 
properties.  As  a  first  step  in  this  direetion,  we  propose  in  this  paper  to  provide  some 
guidelines  for  the  use  of  shock  tube  ignition  data  in  the  development  and  validation  of 
combustion  reaction  mechanisms.  Primarily  we  will  be  coneerned  with  understanding 
the  uneertainties  assoeiated  with  shoek  tube  ignition  data,  as  well  as  the  assumptions  used 
in  modeling  these  data. 

The  use  of  shook  tubes  in  ohemical  kinetio  rate  investigations  has  been  discussed 
before,  e.g.  see  Belford  and  Strehlow  (1969)  [1],  Bowman  and  Hanson  ([2]  and  Lifshitz 
(2001)  [3],  as  well  as  in  detonation  studies,  e.g.  see  Sohultz  and  Shepherd  (2000)  [4]. 

What  is  Ignition  Time? 

Reaction  models  allow  us  to  make  a  direot  comparison  of  the  different  definitions 
of  ignition  time.  Fig.  1  shows  the  model  predictions  for  four  types  of  ignition  time 
indieators  for  five  different  initial  eonditions  that  span  a  wide  region  of  the  initial 
condition  parameter  spaee  for  iso-ootane/02/N2  ignition  mixtures. 

Several  observations  about  the  suitability  of  these  ignition  predictors  can  be  made. 
As  seen  on  a  semi-log  plot,  the  CH*  (excited  state)  emission  has  a  eontinuous  and 
eomplex  behavior,  and  laeks  an  unambiguous  single  feature.  The  CH*  emission  signal 
has  been  simulated  by  assuming  that  it  is  proportional  to  the  produetion  rate  of  CH* 
whieh  is  formed  by  the  reaction  C2H  +  O  ^  CH*  +  CO  and  then  rapidly  quenehed. 
When  measured  using  a  linear  scale,  however,  there  is  usually  a  single  strong  peak  that 
ean  be  identified  with  the  ignition  time.  Pressure  is  a  good  indicator  of  ignition  at  high 
fuel  concentrations  (exeept  in  some  eases  where  the  stoichiometry  is  substantially 
different  than  unity  and  the  ignition  proeess  is  substantially  drawn  out  in  time.)  In  the 
two  low  eoncentration  ignition  cases,  2%  O2,  (|)=1.0  and  0.2%  O2,  (|)=0.25,  however,  there 
is  a  small,  almost  immeasurable  pressure  rise,  whieh  renders  this  ignition  indieator 
unsuitable  for  low  concentration  mixtures.  The  CH*  and  OH  (and  intermediate  species 
CsHe)  mole  fraetion  histories  all  show  elear  evidence  of  a  ehange  owing  to  ignition  for  all 
oases,  yet  there  is  no  unique  way  to  define  ignition  time.  For  example,  the  ignition  delay 
time  could  reasonably  be  defined  as  ocourring  either  at  the  time  of  maximum  rate  of 
ehange  or  of  the  peak  value  of  some  species  or  variable  suoh  as  [OH],  [CH]  or  pressure, 
or  could  be  based  on  an  extrapolation  of  the  maximum  slope  to  the  zero  signal  level. 
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Fig.  1.  CH*  production,  pressure,  OH  and  C3H6  mole  fraction  time  histories  for  a  variety  of  iso- 
octane/02/N2  ignition  mixtures  at  initial  reflected  shock  conditions  of  1400  K,  2  atm:  {A,  20%  O2,  (|)=2.0}, 
{B,  20%  O2,  (|)=0.25},  {C,  0.2%  O2,  (|)=2.0},  {D,  0.2%  O2,  (|)=0.25},  {E,  2%  O2,  (|)=1.0}.  Modeling  using 
the  Davis  and  Law  (1998)  [6]  iso-octane  oxidation  mechanism. 

Table  1  shows  a  comparison  of  the  different  ignition  times  derived  from  several 
different  indicators:  CH*  emission,  pressure  change,  OH  or  C3H6  species  concentration. 
The  variation  of  ignition  times  using  these  definitions,  except  in  case  B,  is  typically  less 
than  2%  from  the  average,  indicating  that  any  of  these  methods  can  be  used  to  identify 
ignition  time  and  that  they  are  reasonably  interchangeable,  over  the  range  of  reaction 
conditions  studied.  In  very  lean  mixtures  such  as  case  B,  20%  O2,  (|)=0.25,  a  more 
specific  definition  may  be  needed,  though  the  scatter  is  still  typically  less  than  10%. 
What  is  especially  important  to  note  is  that  while  there  is  reasonable  agreement  amongst 
the  methods,  some  methods,  such  as  pressure  measurement,  cannot  be  used  over  a  wide 
range  of  conditions  because  of  signal-to-noise  ratio  issues. 
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Case 

XcH*  peak 

XoH  peak 

(dXc3H6/dt)  max 

(dP/dt)  max 

A 

683 

680 

688 

675 

B 

300/380^ 

364 

340 

292 

C 

17330 

17440 

16960 

17200** 

D 

1910 

2020 

1907 

1900** 

E 

4470 

4620 

4640 

4460 

Table  1.  Ignition  times  from  Fig.  1.  ^  first/second  peak;  **  very  small  change  in  P.  All  times  given  in 
microseconds.  Initial  conditions  of  cases  described  in  Fig.  2. 

Several  issues  remain.  First  is  the  issue  of  limited  shoek  tube  test  time.  Though 
methods  are  available  to  extend  reflected  shock  test  time  by  tailoring  driver  gas  mixtures, 
larger  uncertainties  in  temperature  exist  in  shock  tube  measurement  with  test  times 
greater  than  about  2  ms,  and  for  relatively  small  shock  tube  diameters  or  insufficient 
driver  length,  this  upper  test  time  limit  may  be  as  short  as  500  |is  or  less.  Long  measured 
ignition  times,  of  order  several  ms,  should  be  considered  as  lower  limits  of  the  true 
constant  volume  ignition  time,  as  boundary  layer  effects  will  degrade  the  constant 
volume  assumption  by  increasing  the  test  gas  temperature  significantly  at  longer  times. 
This  effect  is  described  in  more  detail  in  a  later  section. 

A  related  issue  is  the  minimum  measurable  ignition  time.  For  combustion 
temperatures,  reflected  shocks  travel  at  a  speed  of  about  0.5  mm/jis,  and  for  optical 
measurements  through  windows  of  the  order  of  6-12  mm  diameter,  this  finite  shock 
transit  time  can  cause  an  experimental  convolution  of  the  signal  of  at  least  12  |a,s  if  the 
viewing  area  is  not  limited  by  the  probe  beam  diameter  or  the  collection  slit  width. 
Piezo-electric  pressure  transducers  typically  have  active  diameters  of  approximately  4 
mm  and  response  frequencies  of  400  kHz,  giving  similar  experimental  convolutions  if 
they  are  mounted  on  the  shock  tube  side  wall,  and  are  prone  to  ringing  if  placed  on  the 
shock  tube  end  wall.  This  ringing  prevents  identification  of  ignition  times  that  are  shorter 
than  several  crystal  ringing  cycles.  Blast  wave  effects  found  in  high  concentration 
mixtures  cause  modifications  to  ignition  times  measured  away  from  the  endwall. 
Corrections  to  ignition  time  measurements  can  be  of  the  order  of  10-40  |is  for  sidewall 
measurement  locations  20  mm  from  the  endwall.  See  Petersen  et  al.  [5]  for  details  of  this 
calculation.  Shock  tube  endwall  emission  measurements  of  ignition  time  can  eliminate 
this  particular  problem.  We  have  found  through  experience  that  measured  ignition  times 
shorter  than  50  |is  have  larger  uncertainties  and  are  less  suitable  for  reliable  study. 

Another  issue  is  multi-stage  ignition.  The  ignition  process  in  large  fuels  such  as 
n-heptane  and  iso-octane  proceed  through  a  series  of  decomposition  and  oxidation  steps 
to  small,  and  more  stable  intermediates,  e.g.  H2,  C2H4,  C2H2,  C3H6,  followed  by  the  final 
oxidation  steps  of  these  small  stable  species.  Example  evidence  of  the  rapid  formation  of 
the  stable  intermediate  C3H6  in  iso-octane  ignition  modeling  is  presented  in  Fig.  1 .  An 
ignition  time  based  solely  on  the  disappearance  of  the  initial  reactant,  in  this  case  iso¬ 
octane,  would  suggest  an  ignition  time  of  the  order  of  10-100  |is,  while  an  ignition  time 
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based  on  the  largest  step  in  pressure,  the  final  removal  of  C3H6,  or  the  rapid  formation  of 
OH  radicals  would  be  two  orders  of  magnitude  longer. 

Depending  on  the  initial  conditions,  the  ignition-indicating  diagnostic  (pressure, 
OH,  etc.)  may  vary  too  slowly  because  of  the  smoothing  effect  of  these  multi-stage 
processes  to  offer  an  obvious  ignition  time.  In  these  cases  an  overall  ignition  time  would 
need  to  be  arbitrarily  defined.  This  ignition  time  should  have  a  kinetic  basis,  for 
example,  time  to  half  the  maximum  change  in  pressure  or  to  half  the  maximum  species 
concentration,  rather  than  a  laboratory  basis,  such  as  time  for  the  CH*  emission  signal  to 
reach  1  volt,  etc.,  which  would  be  difficult  to  relate  to  kinetic  modeling. 

Are  Ignition  Time  Correlations  Useful? 

A  review  of  the  shock  tube  literature  for  a  common  fuel  will  usually  bring  to  light 
a  variety  of  ignition  time  studies.  Shown  in  Fig.  2  are  example  data  from  five  groups 
who  studied  ignition  times  in  n-heptane/Oi/Ar  mixtures.  For  ease  of  comparison,  these 
data  have  been  normalized  using  the  correlation  of  Burcat  et  al.  (1991)  [7].  These  studies 
used  a  variety  of  ignition  time  markers  including:  pressure  rise  and  OH*  and  CH* 
emission.  The  reaction  conditions  vary  over  a  range  of  temperature,  pressure, 
equivalence  ratio  and  fuel  or  oxygen  concentration. 

What  is  clearly  evident  in  the  results  obtained  in  these  different  studies  is  that  they 
can  vary  as  much  as  one  order  of  magnitude  above  and  below  the  correlation  suggested 
Burcat  et  al.  What  is  the  reason  for  this?  Is  this  caused  by  use  of  a  correlation  that  is  not 
valid  for  the  conditions  of  the  other  studies?  (As  will  be  evident  from  the  modeling 
discussion  below,  current  correlations  can  still  give  a  fair  first  approximation  to  the 
ignition  times,  even  outside  their  data  set.)  Is  this  caused  by  differences  in  the  definition 
of  ignition  time  amongst  the  studies?  (Apparently  not,  based  on  the  example  (Fig.  1, 
Table  1)  discussed  above.)  Is  this  caused  by  errors  in  the  experimental  methods  of  the 
different  workers?  (This  type  of  error  is  really  only  detectable  through  comparisons  based 
on  correlations.)  Is  this  caused  by  differences  in  the  shock  tubes?  (The  performance  of 
particular  shock  tubes,  which  is  often  related  to  their  diameter,  may  not  be  suitable  for 
ignition  time  studies.)  Answers  to  these  questions  are  needed  before  these  data  can  be 
used  for  quantitative  comparison  with  computational  models.  Resolution  to  these 
questions  will  hopefully  enable  the  modeler  to  use  shock  tube  ignition  time  data  to 
constrain  the  predictions  of  their  reaction  mechanism  to  better  than  the  nearly  two  orders 
of  magnitude  variation  evident  in  this  plot. 

The  development  of  correlations  for  ignition  time  data  is  necessary  for  several 
reasons.  Firstly,  some  modelers  require  ignition  time  correlations  to  embed  in  their 
engine  programs  or  reactive  flow  combustion  codes  rather  than  use  full  or  reduced  kinetic 
modeling  to  describe  the  chemical  processes.  Correlations  also  guide  the  experimentalist 
in  the  design  of  data  sets,  by  reducing  the  number  of  experimental  conditions  needed  to 
fully  examine  the  ignition  behavior  of  a  particular  fuel.  Additionally,  ignition  time 
correlations  are  essential  to  allow  comparison  of  findings  from  different  studies  that 
might  have  been  conducted  with  different  reaction  conditions.  Such  comparisons  are 
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necessary  as  can  be  seen  from  Fig.  2,  because  some  data  sets  fall  far  from  the  consensus 
and  may  not  be  consistent  with  other  data  sets,  and  may  yield  misleading  results  if  used 
in  model  validation.  Comparisons  such  Fig.  2  also  typically  provide  a  means  of  assessing 
the  degree  of  data  scatter  in  individual  studies. 


Fig.  2.  n-Heptane  ignition  times.  All  times  correlated  using  scaling  of  Burcat  et  al.  (1991):  Tign  varying  as 
[C7H16]  ■’■'[Ar]°^ 


Correlations  can  take  many  forms.  Over  the  range  of  conditions  that  Burcat  et  al. 
measured  n-heptane  ignition  times,  the  authors  correlated  ignition  times  using  an 
expression  of  the  form; 

x,g„  =  C  [C7H16]'’  [02]"^  [Ar]’^  exp(EAB/RT)  Eqn.  1 . 

This  form,  an  improvement  over  earlier  expressions  without  a  dependence  on  bulk  carrier 
gas  concentration,  describes  the  behavior  of  n-heptane  ignition  times  using  five 
parameters,  and  assumes  that  the  ignition  times  vary  as  power  laws  in  the  concentrations 
of  fuel,  oxygen  and  argon,  respectively,  and  exponentially  as  (1/T).  See  Eifshitz  (2001) 
[3]  for  a  discussion  about  this  form  of  correlation. 

Recent  work  in  our  laboratory  by  Horning  et  al.  (2002)  [8]  found  that  for  n- 
alkanes,  ignition  times  could  be  correlated  well  with  an  expression,  also  with  five 
parameters,  of  the  form: 

Tign  =  C  Xo2“  exp(EAH/RT)  Eqn.  2. 

and  in  particular  for  n-heptane; 
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Xign  =  6.67x10'^^  exp(44,600/RT)  Eqn.  3. 

where  x  ign  is  in  s,  T  is  in  K,  P  is  in  atm,  and  Ea  is  in  eal/mol. 

The  key  advantage  of  this  form  of  the  correlation  is  that  it  mimics  several  current 
n-heptane  oxidation  mechanisms  in  the  dependence  of  ignition  time  on  P,  X02  and  Ea  at 
temperatures  above  1100  K.  Over  the  limited  equivalence  ratio  range,  (|)  =  0.5  to  2.0,  a 
power  law  relation  provides  a  good  simple  description  of  the  dependence  of  ignition  time 
with  equivalence  ratio.  Eor  a  wider  equivalence  ratio  range,  (|)  =  0.2  to  4.0,  a  different 
form  of  the  correlation  is  preferred,  one  that  duplicates  the  experimental  and  model 
predictions  of  a  weak  power  law  dependence  at  small  (|),  where  oxygen  chemistry  plays  a 
more  important  role  in  determining  ignition  time,  and  a  strong  power  law  dependence  at 
large  (|),  where  fuel  pyrolysis  and  oxygen  scavenging  chemistry  plays  an  important  role. 
The  Horning  et  al.  correlation  takes  advantage  of  the  predictions  and  knowledge  derived 
from  a  fuller  kinetic  model  of  the  ignition  process.  (See  Horning  et  al.  (2002)  for  a 
comparison  of  the  ignition  time  predictions  of  several  models  and  their  relationship  with 
this  correlation.) 

Eor  the  case  of  iso-octane,  Davidson  et  al.  (2002)  [9]  found  ignition  time 
dependences  that  were  similar  to  those  Homing  observed  in  n-heptane  for  pressure, 
equivalence  ratio,  and  activation  energy,  but  found  a  stronger  dependence  on  fuel  or 
oxygen  mole  fraction  at  large  concentrations.  The  correlation  proposed  by  Davidson  et 
al.  for  the  high  concentration  data  is 

Xign  =  4.50x10'^°  exp(-232  Xpuei)  exp(44,780/RT).  Eqn.  4. 

Another  form  of  the  correlation  that  captures  the  very  low  concentration  iso-octane  data, 
gives  similar  high  concentration  predictions,  and  is  related  to  the  Homing  et  al. 
correlation  is  given  by: 

Xign  =  1.65x10'^°  Xo2'°-^^  exp(43,050/RT)  exp(-129  Xm).  Eqn.  5. 

These  two  examples,  the  non-power  law  variation  of  n-heptane  ignition  time  with 
equivalence  ratio  and  the  non-power  law  variation  of  iso-octane  ignition  time  with  fuel 
concentration,  demonstrate  that  a  simple  power  law  correlation  for  ignition  time 
measured  over  a  very  wide  range  of  conditions  may  not  always  capture  the  needed 
chemistry  and  phenomena.  Correlations  using  power  laws  are  widespread  in  the 
literature,  and  for  data  sets  of  limited  extent  in  pressure,  equivalence  ratio,  fuel 
concentration  and  temperature,  may  still  be  a  very  good  choice.  A  comparison  with 
detailed  kinetic  model  predictions  is  usually  warranted  to  understand  the  useable  regime 
of  conditions  for  a  correlation. 
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Comparing  Correlations  and  Kinetic  Modeling 

As  discussed  in  the  previous  section,  it  is  critical  when  eomparing  ignition  time 
data  with  eorrelations  and  kinetic  models  to  know  how  eomplete  the  experimental 
eoverage  of  the  (T,  P,  X02,  (|))  parameter  spaee  is  for  a  particular  fuel  and  buffer  gas: 
M=Ar  or  N2.  In  the  following  discussion  we  will  eontinue  to  use  as  examples  the  n- 
heptane  ignition  times  by  Horning  et  al.  (2002)  and  iso-oetane  ignition  times  by  Davidson 
et  al.  (2002). 

Homing  et  al.  measured  n-heptane  ignition  times  using  endwall  CH*  emission 
over  a  range  of  eonditions  that  ineluded  1329-1676  K,  1.15-5.69  atm,  0.2-1.75  %  n- 
heptane,  2.2-19.7  %  oxygen  and  equivalenee  ratios  between  0.5  and  2.0.  The  measured 
ignition  times  for  these  experiments  fell  between  100  and  500  |is  and  minimized  any 
error  assoeiated  with  boundary-layer-generated  temperature  non-uniformity  in  the 
relatively  large  diameter  (15.3  em)  shoek  tube.  Similarly  Davidson  et  al.  measured  iso- 
oetane  ignition  times  over  a  range  of  eonditions  that  ineluded  1177-2009  K,  1.18-8.17 
atm,  lOOppm  -  1%  iso-oetane,  0.125  to  12.5%  oxygen  and  equivalence  ratios  of  0.25  to 
2.0.  Work  by  Niemitz  et  al.  (1988)  [10]  extends  this  range  to  1.6%  iso-oetane  and  20% 
oxygen. 

The  extent  of  parameter  space  that  these  data  eover  in  (|)-Xo2  spaee  is  indieated  in 
Fig.  3.  In  both  oases,  eoverage  of  the  (|)-Xo2  spaee  is  good,  but  it  is  evident  when 
presented  in  this  manner  that  further  measurements  at  high  oxygen  oonoentrations  (20%) 
and  equivalenee  ratios  of  0.25  and  2.0  may  be  warranted.  The  extent  of  the  eoverage  in 
pressure  and  temperature  space  is  less  important  for  these  eonditions  (1100-2000  K  and 
1-8  atm)  beoause  both  the  experimental  data  and  ourrent  models  show  simple  power  law 
variation  of  the  ignition  times  with  pressure  and  exponential  variation  with  inverse 
temperature,  and  hence  extrapolation  to  other  conditions  can  be  made  more  eonfidently. 

We  ean  eompare  the  experimental  ignition  time  data  (in  the  form  of  a  eorrelation) 
with  the  predietions  of  a  kinetie  model  for  a  wide  range  of  eonditions  by  using  a  three- 
dimensional  representation  over  the  (|)-Xo2  spaee.  In  the  top  frame  of  Fig.  4  the  ignition 
time  predietions  from  the  Davidson  et  al.  iso-oetane  shoek  tube  eorrelation  are  shown 
along  with  the  eonstant  volume  model  predietions  using  the  Pitseh  et  al.  (2001)  [11] 
redueed  reaetion  mechanism.  This  three-dimensional  plot  shows  the  variation  of  ignition 
time  with  oxygen  mole  fraetion  (from  100  ppm  to  20%)  and  equivalenee  ratio  (from  0.25 
to  2.0)  for  initial  eonditions  of  1400  K  and  2  atm.  The  bottom  frame  shows  similar 
results  from  the  Horning  et  al.  n-heptane  data  and  the  Pitseh  et  al.  (2000)  [12]  redueed 
reaetion  mechanism.  The  advantage  of  this  representation  is  that  a  direet  eomparison  of 
eorrelation  (and  henee  data)  and  model  ean  be  made  over  wide  range  of  eonditions. 
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Fig.  3.  Extent  of  parameter  space  covered  by  experimental  data  of  Homing  et  al.  for  n-heptane  ignition 
time  measurements  and  Davidson  et  al.  iso-octane  data.  Open  circle  in  iso-octane  chart  is  data  of  Niemitz 
etal.  (1988). 

We  see  that  the  Pitseh  et  al.  iso-oetane  model  (47  species,  134  reactions)  captures 
the  major  features  of  the  iso-octane  data  correlation:  the  variation  with  stoichiometry  and 
the  strong  variation  with  oxygen  mole  fraction.  This  particular  correlation  suggests  that 
at  large  oxygen  mole  fraction  20%  and  large  equivalence  ratio  2.0  that  the  ignition  time 
shortens  dramatically.  This  region  of  mixture  space  in  Fig.  3  has  not  been  tested 
experimentally,  but  the  Pitseh  et  al.  model  does  suggest  that  it  will  still  be  shortened 
somewhat. 

In  the  case  of  n-heptane,  the  ignition  time  correlation  is  very  simple  (power  law  in 
both  oxygen  mole  fraction  and  equivalence  ratio).  But  the  kinetic  model  of  Pitseh  et  ah, 
though  able  to  generate  approximately  correct  ignition  times,  does  not  show  the  same 
simple  trends  over  the  whole  surface.  Anomalous  non-monotonic  variation  with  oxygen 
concentration  at  lean  equivalence  ratios  is  evident,  for  example.  The  present 
experimental  database  for  lean  mixtures  is  currently  not  extensive  enough  to  support  or 
contradict  the  variation  found  in  the  model.  This  approximate  representation  of  the 
ignition  time  correlation  (and  hence  data)  is  typical  of  a  strongly  reduced  reaction 
mechanism,  in  this  case,  44  species  and  112  reactions.  At  large  oxygen  mole  fraction 
20%  and  small  equivalence  ratio  0.25,  the  correlation  and  the  model  diverge 
significantly.  This  mixture  regime  has  also  not  been  tested  experimentally  based  on  the 
n-heptane  parameter  space  coverage  of  Fig.  3. 

One  observation  directly  evident  from  the  above  presentation  format  is  that  at  a 
high  oxygen  mole  fraction  of  20%  and  equivalence  ratios  of  approximately  1,  the 
correlations  and  models  predict  that  the  iso-octane  ignition  times  will  be  shorter  than  the 
n-heptane  ignitions.  This  trend  is  not  expected  based  on  assumed  role  of  iso-octane  in 
reducing  knock  (i.e.  extending  ignition  times)  in  gasoline  engines. 

This  observation  that  the  iso-octane  and  n-heptane  ignition  times  are  similar  is  not 
new.  In  Westbrook,  Warnatz  and  Pitz  (1988)  [10]  there  is  evidence  (their  Fig.  I)  that  the 
ignition  time  of  iso-octane  is  comparable  with  that  of  n-heptane  under  certain  high 
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concentration  conditions.  Westbrook  et  al.  (2002)  [13]  have  also  shown  in  rapid 
eompression  maehines  that  it  is  possible  to  have  non-monatonie  behavior  in  the 
relationship  between  oetane  number  and  knoeking,  higher  oetane  number  fuels  produeing 
more  knoeking  than  lower  oetane  number  fuels  in  some  eases.  The  differenees  in  the 
ehemieal  kinetie  history  experieneed  in  shook  tubes  and  rapid  eompression  maehines  and 
in  praotioal  engines  is  also  disoussed  in  a  review  by  Westbrook  et  al.  (2000)  [14]. 


Fig.  4.  Three-dimension  comparison  of  correlations  and  models  for  n-heptane  and  iso-octane  ignition 
times. 


An  extensive  shook  tube  ignition  time  study  oan  offer  more  than  just  a  series  of 
disorete  targets  for  kinetie  modelers.  Multi-dimensional  representations  of  the  data  oan 
expose  global  trends  which  themselves  oan  be  used  as  kinetie  targets  for  model  validation 
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and  refinement.  Correlations  offer  a  very  simple  way  to  do  this,  making  signifieant 
similarities  and  differences  and  important  trends  in  the  experiment  and  the  model  more 
evident. 

Shock  Tube  Constant  Volume  Behavior 

In  comparing  modeled  ignition  data  with  experimental  data,  the  model 
calculations  are  usually  done  with  a  constant  volume  or  constant  density  constraint.  One 
measure  of  the  validity  or  suitability  of  this  assumption  is  a  comparison  of  modeled  and 
experimental  pressure  jump  and  plateau  pressures  during  and  after  ignition. 

Example  shock  tube  pressure  traces  for  a  range  of  fuel  concentrations  are  shown 
in  Fig.  5.  These  traces  have  several  identifiable  features.  The  first  plateau  is  the  filling 
pressure  Pi,  the  second  plateau  is  the  pressure  behind  the  incident  shock  P2.  The  pressure 
transducer  in  these  experiments  was  placed  2  cm  from  the  end  wall  of  the  shock  tube, 
resulting  in  a  short  interval,  of  order  65  |is  where  this  incident  shock  pressure  can  be 
measured.  The  third  plateau  is  the  reflected  shock  pressure  P5.  The  rapid  rise  in  pressure 
at  ignition  in  the  highest  concentration  examples,  0.5%  and  0.25%  fuel,  seen  in  both  the 
model  and  data  of  the  upper  frames  of  Fig.  5,  is  the  constant  volume  ignition  pressure 
increase.  The  ragged  form  of  the  experimental  pressure  data  after  this  time  is  evidence  of 
a  blast  wave  or  detonation  that  occurs  after  this  pressure  jump.  In  the  lower 
concentration  examples,  the  post-ignition  plateau  pressures  are  similar  to  P5  and  the 
ignition  pressure  jump  is  barely  noticeable.  At  times  longer  than  2  ms  in  all  of  these 
traces,  the  pressure  begins  to  fall  because  of  the  interaction  of  the  driver  rarefaction  wave 
with  the  reflected  shock  region. 

If  the  shock  tube  were  an  ideal  constant  volume  reactor,  then  even  with  the  energy 
release  process  occurring  during  ignition,  the  modeled  pressure  and  measured  pressure 
would  achieve  the  same  levels  during  the  initial  ignition  pressure  jump  and  the  final 
resting  plateau.  The  ignition  pressure  jump  is  determined  by  the  energy  release  rate  of 
the  ignition  process,  and  the  final  resting  pressure  is  determined  by  the  energy  that  is 
released.  A  comparison  of  the  experimental  and  modeled  ignition  pressure  jumps  is 
shown  in  Fig.  6.  At  higher  fuel  concentrations,  the  pressure  jump  is  substantially  less 
than  that  predicted  by  the  constant  volume  model,  which  is  indicative  of  a  failure  of  the 
constant  volume  constraint  during  the  last  stages  of  these  energetic  ignition  events. 
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Fig.  5.  Side  wall  pressure  measurements  and  modeling  during  ignition  events  in  iso-octane/02/Ar  for  a 
variety  of  fuel  concentrations;  0.5%  to  0.01%.  Initial  shock  conditions:  ~1500K,  -1.4  atm,  (|)=1.0. 
Modeling  using  Davis  and  Law  (1998). 


Fig.  6.  A  comparison  of  the  measured  and  modeled  pressure  jump  ratio  and  plateau  pressure  ratio  during 
ignition  events  in  iso-octane/02/ Ar,  ~1500K,  1. 4-5.0  atm,  (|)=1.0. 


Because  this  initial  pressure  jump  occurs  exponentially,  its  influence  is  limited  to 
a  small  fraction  of  the  time  prior  to  ignition.  Hence,  the  measured  pressure  profile  of  the 
ignition  process  can  be  approximated  by  the  calculated  pressure  profde  using  a  constant 
volume  constraint  up  to  the  time  of  ignition.  After  that  time,  for  ignition  events  with 
large  energy  release,  the  shock  tube  does  not  behave  like  a  constant  volume  reactor. 
Figure  6  shows  that  only  for  stoichiometric  iso-octane  ignition  experiments  with  mixtures 
of  '-'500  ppm  or  less  of  fuel  can  we  assume  that  the  entire  ignition  process,  including  the 
plateau  region  after  the  ignition,  can  be  accurately  modeled  using  a  constant  volume 
constraint. 

One  application  of  this  pressure  jump  analysis  is  determining  which  shock  tube 
ignition  experiments  can  be  quantitatively  modeled  throughout  the  entire  ignition 
process.  Thus  in  the  case  of  0.05%  iso-octane,  we  may  confidently  use  the  constant 
volume  constraint  to  model  quantitative  species  concentration  measurements  made  after 
the  point  of  ignition.  An  example  of  OH  data  from  such  a  shock  tube  experiment, 
Davidson  et  al.  (2002)  [9],  is  shown  in  Fig.  7.  These  OH  data  are  based  on  cw  ring-dye 
laser  absorption  measurements  of  the  OH  A-X  (0,0)  band  Ri(7)  line;  concentration 
uncertainties  are  of  the  order  of  ±5%. 

In  this  experiment  the  difference  between  the  final  plateau  values  of  the  model 
and  data  can  be  attributed  to  discrepancies  in  the  model.  Either  inaccuracies  in  the 
reaction  rate  coefficients  or  thermo-chemical  data  could  be  responsible.  Recent  work  by 
Herbon  et  al.  (2003)  [15]  has  indicated  that  some  of  this  variation  in  combustion  product 
partial  equilibrium  values  can  be  attributed  to  errors  in  the  value  used  for  the  OH  heat  of 
formation. 


Fig.  7.  Side  wall  OH  absorption  measurements  and  modeling  during  ignition  events  in  iso-octane/02/Ar, 
Davidson  et  al.  (2002)  [9]. 


Boundary  Layer  Effects  on  Shock  Tube  Temperatures 

The  following  discussion  is  based  on  corrections  to  shock  tube  ignition  times 
given  by  the  boundary  layer  model  of  Petersen  and  Hanson  (2002)  [5]. 
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In  the  ideal  deseription  of  a  shoek  wave  experiment,  the  ineident  shoek  front 
passes  along  the  length  of  the  shoek  tube,  leaving  behind  a  region  of  shoeked  gas  with 
uniform  P2,  T2,  and  V2,  the  gas  veloeity  in  laboratory  eoordinates.  This  is  followed  by  a 
refleetion  of  the  shoek  wave  from  the  end  wall  of  the  shoek  tube  and  formation  of  a 
uniform  refleeted  shook  region  with  uniform  P5,  T5  and  V5=0. 

The  actual  behavior  of  these  shock  waves  is  modified  by  formation  of  boundary 
layers  behind  the  incident  and  reflected  shocks.  The  result  of  the  boundary  layer  behind 
the  incident  shock  is  to  generate  a  gradul  increase  in  temperature  and  pressure  in  the  core 
of  the  flow.  When  the  reflected  shock  passes  through  this  non-uniform  pressure  and 
temperature  distribution,  it  also  experiences  an  increase  in  temperature  and  pressure,  and 
this  increase  appears  to  be  linear  with  test  time. 

The  effects  of  boundary  layers  on  shock  tube  performance  are  unavoidable, 
though  it  is  possible  to  reduce  these  effects  to  a  negligible  level  in  some  cases.  The 
primary  ways  to  reduce  the  influence  of  boundary  layers  are  to  employ  a  large  shock  tube 
diameter,  raise  the  test  gas  pressure  or  decrease  the  portion  of  the  test  time  that  is  used  in 
the  measurement.  One  simple  way  to  quantify  this  boundary  layer  effect  is  by 
considering  the  value  dT*/dt  which  equals  the  relative  change  of  T*  =  AT5/T5  initial  with 
reflected  shock  test  time.  Typically,  large  diameter  shock  tube  have  dT*/dt  values  less 
than  10  s'\  This  effect  will  also  be  a  function  of  the  initial  fill  pressure  Pi.  Thus  a  6-inch 
diameter  (15  cm)  shock  tube  could  have  a  change  of  temperature  attributed  to  boundary 
layer  effects  over  1  ms  of  test  time  of  less  than  15  K  for  a  1500  K  reflected  shock.  For  a 
shock  tube  less  than  4  inches  diameter  (10  cm)  this  effect  can  rapidly  increase  to  values 
of  20  to  30  s'^  or  more. 

The  increase  in  temperature  due  to  the  boundary  layer  has  a  dramatic  effect  on 
apparent  activation  energy  of  ignition  times.  Fig.  8  and  Table  2  illustrate  the  effect  that  a 
shock  tube  with  a  small  diameter  (and  hence  a  large  dT*/dt)  and  an  experimental  data  set 
with  long  (1-2  ms)  test  times  have  on  a  hypothetical  ignition  time  measurement  where  the 
true  activation  energy  is  40  kcaFmol.  Substantial  errors  in  the  apparent  activation  energy 
(greater  than  6  kcal/mol)  are  possible. 


16 


Fig.  8.  Simulated  ignition  times  showing  the  effect  of  boundary-layer-induced  test  gas  heating  on  the 
ignition  time  activation  energy  for  a  hypothetical  fuel  with  an  ignition  time  activation  energy  of  40 
kcal/mol. 

Without  a  quantification  of  the  value  of  dT*/dt  for  a  small  diameter  shoek  tube 
(an  effeet  which  is  pressure  and  temperature  dependent),  eorrections  to  previously 
published  ignition  time  measurements  are  diffieult.  However,  an  understanding  of  the 
signifieanee  of  dT*/dt  ean  be  useful  in  explaining  the  apparent  differenees  in  aetivation 
energy  found  in  shook  tubes  of  different  diameters. 


Shook  Tube 
Performanoe 

Inoluding  ignition 
times  up  to  500  |is 

Inoluding  ignition 
times  up  to  1  ms 

Inoluding  ignition 
times  up  to  2  ms 

dT*/dt 

[s-'] 

Measured  Ea 
[oal/mol] 

Measured  Ea 
[oal/mol] 

Measured  Ea 
[oal/mol] 

0 

40000 

40000 

40000 

10 

39200 

38400 

37600 

20 

38600 

37200 

35800 

30 

37800 

35600 

33500 

Table  2.  Example  effect  on  ignition  activation  energy  by  boundary-layer-induced  shocked  gas  heating. 
Measured  Ea  for  true  ignition  activation  energy  of  40  kcakmol  and  data  over  the  range  of  1 176  to  1428K. 

Effect  of  Bulk  Carrier  Gas 

To  improve  shook  tube  performanoe,  shook  wave  experiments  are  normally 
performed  in  argon,  rather  than  nitrogen.  This  eliminates  the  possible  infiuenoe  of 
vibrational  relaxation  in  the  oarrier  gas,  as  well  as  eliminating  or  minimizing  shook  wave 
bifuroation  near  the  wall  in  the  refleoted  shook  front.  As  well,  for  the  same  initial  shook 
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pressure  Pi,  argon  shocks  can  have  higher  reflected  shock  temperatures,  as  no  energy  has 
to  go  into  the  vibrational  modes  of  the  shock-heated  gas.  Argon  shocks  have  historically 
been  considered  to  have  “quieter”  pressure  profiles  and  this  may  be  related  to  differences 
in  boundary  layer  growth. 

Normal  combustion  occurs  in  air,  of  course,  and  differences  in  ignition  time  data 
collected  in  argon  and  air  should  be  recognized  and  quantified,  as  should  differences  in 
modeling  with  these  two  different  carrier  gases.  Figure  9  shows  modeling  of  iso-octane 
ignition  using  two  different  reduced  models  for  initial  conditions  of  1400  K,  2  atm,  2% 
O2,  (|)  =  1.0  with  argon  and  nitrogen  as  the  carrier  gas.  For  one  model,  Davis  and  Law, 
the  differences  are  insignificant,  while  for  the  other  model  by  Pitsch  et  al.  the  difference 
is  significant. 

A  comparison  between  ignition  times  in  gas  mixtures  with  nitrogen  or  argon  as 
carrier  gases  allows  the  investigator  to  study  the  effect  of  small  temperature  changes 
during  the  ignition  process.  During  the  initial  fuel  decomposition  step,  the  temperature  of 
the  shock  heated  gas  mixture  drops  slightly  as  energy  is  extracted  to  break  fuel  bonds. 
This  temperature  drop  is  larger  in  argon  than  in  nitrogen.  The  Pitsch  et  al.  mechanism 
shows  much  more  rapid  growth  of  the  temperature  in  the  argon  carrier  gas,  even  though 
the  argon  carrier  gas  minimum  temperature  is  lower  than  that  found  with  the  nitrogen 
carrier  gas.  Once  the  temperature  of  the  secondary  reaction  products  begins  to  increase  it 
rapidly  accelerates  the  ignition  process.  In  the  Davis  and  Law  mechanism  this  increase  is 
not  evident  in  the  M=Ar  case,  indicating  that  the  sub-mechanism  controlling  the 
secondary  reaction  products  is  not  as  strongly  temperature  sensitive.  Evidence  exists  for 
some  fuels  that  there  is  no  experimental  sensitivity  to  carrier  gas  but  not  for  all  fuels.  See 
Petersen  et  al.  (1999)  [16]  for  a  discussion  of  this  effect  in  methane  oxidation. 

The  effect  of  the  heat  capacity  of  the  bulk  carrier  gas  on  shock  tube  test  gas 
temperature  is  different  than,  but  related  to  the  important  effect  of  the  time  scale  of 
vibrational  relaxation  of  a  diatomic  carrier  gas  (i.e.  nitrogen  and  oxygen)  on  the  shock 
tube  test  gas  temperature.  For  very  dilute  mixtures  ('-'0.1  %  fuel)  in  nitrogen  or  air  at  low 
temperatures  and  pressures,  vibrational  relaxation  times  can  be  substantially  longer  than 
the  reflected  shock  test  time.  In  these  cases,  vibrationally-frozen  thermodynamic 
properties  should  be  used  for  determining  shock  temperatures  and  pressures.  At  high  fuel 
concentrations,  high  pressures  and/or  high  temperatures,  vibrational  relaxation  times  are 
usually  short  enough  to  use  equilibrium  thermodynamic  properties  in  the  shock  equations 
for  determining  shock  temperatures  and  pressures.  Simple  solutions  for  the  intermediate 
regime  have  not  been  developed,  and  a  standard  treatment  of  this  combination  or 
coupling  of  the  vibrational  relaxation  equations  with  the  kinetic  equations  is  needed. 
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Fig.  9.  A  comparison  of  modeled  iso-octane  ignition  using  two  different  carrier  gases.  Initial  conditions  of 
1400  K,  2  atm,  2%  O2,  (|)  =  1.0.  Left  frame,  Davis  and  Law  mechanism  [6];  right  frame  Pitsch  et  al. 
mechanism  [11]. 

Effect  of  Impurities  on  Measurements 

Finally,  it  should  be  mentioned  that  impurities  in  the  initial  mixture  or  shook  tube 
itself,  oan  in  some  oases  have  a  signifioant  effeot  on  the  ignition  time.  Table  3  shows 
predicted  ignition  times  from  model  calculations  for  a  variety  of  hydrocarbons  under 
typical  shock  tube  conditions.  Impurities  are  simulated  here  by  the  addition  of  1  ppm  H 
atoms  to  the  initial  mixture.  What  is  immediately  evident  is  that  the  effect  of  impurities 
on  various  fuel  species  fall  into  two  groups. 

The  ignition  times  of  ethylene,  acetylene  and  hydrogen  are  strongly  perturbed  by 
the  addition  of  small  amounts  of  H-atoms  to  the  initial  mixtures.  For  these  fuels,  the 
ignition  time  is  dominated  by  the  time  it  takes  for  the  exponential  growth  in  the  radical 
pool.  Increasing  the  H-atom  population,  shortens  this  time. 

In  the  other  fuels,  the  ignition  time  can  be  seen  as  a  multi-step  process,  with  the 
first  step  being  the  rapid  decomposition  of  the  fuel,  and  the  second  step  being  the  slow 
decomposition  of  the  intermediate  products  and  exponential  growth  of  the  radical  pool. 
The  initially-added  H-atoms  are  rapidly  scavenged  in  the  first  step  by  the  fuel  or  the 
major  fuel  decomposition  products  to  form  less  reactive  intermediate  products,  or  by 
reaction  with  O2  that  forms  OH,  which  is  also  rapidly  scavenged  by  the  intermediate  fuel 
products.  The  overall  ignition  time  is  thus  not  strongly  affected. 

It  is  expected  from  this  analysis  that  ignition  time  measurements  of  ethylene, 
acetylene  and  hydrogen  should  show  large  scatter  unless  great  care  is  taken  in 
maintaining  high  purity  conditions  in  the  shock  tube.  The  scatter  in  ignition  times 
measurements  of  larger  hydrocarbons,  on  the  other  hand,  should  not  be  dominated  by  the 
effects  of  impurities. 
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Fuel 

T5  [K] 

Tign  with 
no  H 
added 
[ns] 

Tign  with 
1ppm  H 
added 
[ns] 

% 

change 

Mechanism 

C2H4 

1300 

429 

274 

36.1 

Laskin  et  al.  (2000)  [17] 

C2H2 

1150 

467 

333 

28.7 

Laskin  et  al.  (2000) 

H2 

1050 

244 

191 

21.7 

Smith  et  al.  (1999)  [181 

CD 

X 

0 

1300 

716 

699 

2.4 

Laskin  et  al.  (2000) 

C2H6 

1350 

445 

440 

1.1 

Smith  et  al.  (1999) 

CH4 

1800 

314 

313 

0.3 

Smith  et  al.  (1999) 

C3H8 

1500 

555 

554 

0.1 

Smith  et  al.  (1999) 

n-C/H-ie 

1500 

269 

269 

0.1 

Held  etal.  (1997)  [19] 

iso-CsH-is 

1500 

684 

684 

0.0 

Davis  and  Law  (1998)  [6] 

Table  3:  Effect  of  impurities  on  shock  ignition  times.  All  examples  are  for  5%  O2,  (|)=1,  1  atm,  and 
constant  volume  calculation.  Temperatures  where  selected  so  that  the  modeled  ignition  time 
[^{XoHl/ti^tlmax  wcrc  approximately  500  ps. 

Future  Challenges  and  Needs  for  Shock  Tube  Research 

In  reviewing  the  status  of  ignition  time  measurements  in  shoek  tubes,  several 
ehallenges  and  needs  are  apparent.  The  needs  are  those  of  the  kineties  modeling 
eommunity,  whieh  seeks  to  build  ever  more  aeeurate  and  eomplete  detailed  models  of 
eombustion  ehemistry,  ineluding  ignition  times.  The  key  ehallenges  for  the 
experimentalist  are  to  improve  the  quality  and  extent  of  ignition  data  (e.g.,  through 
measurement  of  relevant  speeies  time  histories  over  a  broader  range  of  mixtures  and  test 
eonditions),  and  to  find  means  to  extend  the  useful  operating  envelope  of  shook  tube 
experiments  outside  the  limited  regime  where  the  ZOAM  is  applioable. 

As  a  speoifio  example,  there  oontinues  to  be  a  need  to  measure  longer  ignition 
times,  those  of  order  1-2  ms  and  longer.  For  these  measurements  to  be  useful,  oonfidenoe 
in  the  temperature  profile  of  the  refieoted-shook-heated  test  gas  is  needed.  Several 
approaohes  are  possible.  One  method  would  be  to  measure  dT(t)/dt  (i.e.  T5(t)  for  eaoh 
experiment.  This  oharaoterization  oan  be  done  using  infrared  emission  measurement 
sohemes  for  non-reaotive  gas  mixtures,  but  this  method  may  not  be  suitable  for  ignition 
experiments.  A  more  praetieal  seheme  would  be  to  relate  refieeted  shoek  temperature 
variation  with  pressure  variation  and  take  advantage  of  the  existing  eapability  to  measure 
pressure  time  histories  during  shoek  tube  experiments.  This  approaeh  ean  be  rigorously 
developed  onee  aeeurate  time-varying  temperature  measurement  teehniques  are 
established  for  reaetive  flows.  However,  temperature  measurements  in  reaetive  flows  in 
shook  tubes  are  diffioult;  see  Sohulz  et  al.  (2002)  [20]  for  one  attempt  to  do  this  with  UV 
CO2  absorption. 

Current  modeling  work  of  engine  eombustion,  in  partioular  homogeneous  oharge 
oompression  ignition  (HCCI)  studies,  requires  ignition  times  for  low  temperature  (700- 
1000  K),  high  pressure  (20-60  atm),  and  rieh  mixtures  (high  fuel  eoneentration  1-2%)  of 
praetieal  fuels  (gasoline),  surrogate  fuels  (binary  and  ternary  mixtures)  and  single  fuel 
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components  (n-alkanes,  branched  alkanes,  cyclo-alkanes,  alkenes,  and  aromatics).  Very 
little  data  of  this  type  exists,  though  early  work  by  Ciezki  and  Adomeit  (1993)  [21]  and 
Fieweger  et  al.  (1997)  [22]  does  start  to  address  some  of  these  needs.  These 
measurements  are  difficult;  in  particular,  measurements  near  700  K,  and  20  atm  are 
expected  to  have  very  ignition  times,  so  that  methods  to  confirm  the  temperature-time 
profde  are  also  needed.  Tailored  driver  gas  mixtures  can  also  be  used  to  extend  shock 
tube  test  time,  but  again  temperature-time  profiles  for  the  potentially  very  long  test  times 
must  be  verified. 

The  use  of  shock  tubes  for  acquiring  species  concentration  time-history  data 
during  energetic  ignition  experiments  is  limited  at  present  in  two  ways.  Either  the  full 
ignition  time-history  can  be  studied  as  a  constant  volume  process  if  the  energy  release  is 
low  (as  in  the  iso-octane  example  described  previously),  or  only  the  early  time  regime 
(before  the  final  large  energy  release)  can  be  studied  if  the  constant  volume  constraint  is 
compromised  during  ignition.  Improved  modeling  of  the  flow  fields  in  shock  tubes 
during  energetic  shocks,  the  development  of  one-dimensional  temperature,  pressure  and 
velocity  profiles  behind  reflected  shocks  for  example,  would  permit  extension  of  the 
range  of  conditions  that  concentration  time-history  data  could  be  usually  acquired  and 
analyzed.  Efforts  are  currently  being  made  to  develop  these  fluid  codes  with  complete 
chemistry  to  permit  a  more  rigorous  validation  of  kinetic  mechanisms  using  species 
concentration  time-history  targets  over  a  wider  regime. 

And  finally,  critical  review  of  existing  and  future  ignition  time  measurements  are 
needed.  Eor  the  reasons  described  earlier,  not  all  published  data  sets  will  necessarily  be 
found  to  be  consistent  with  one  another.  Before  ignition  time  measurements  should  be 
used  in  validating  or  refining  kinetic  mechanisms,  some  estimate  of  the  accuracy  or 
reliability  of  the  data  set  should  be  made.  Reliance  on  multiple  data  sets  from  several 
different  laboratories  for  tuning  reaction  mechanisms  will  enable  modelers  to  place  a 
great  deal  of  scientific  confidence  in  what  continues  to  be  an  important  and  needed 
experimental  activity. 
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